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Abstract 

The hadronization structure of t~ —>■ rj-K^-K°v T decays is analyzed 
within the framework of Resonance Chiral Theory, considering only the 
contribution of the lightest meson resonances at leading order in the 1 /Nc 
expansion. After imposing the asymptotic behavior of vector spectral 
functions ruled by QCD, unknown effective couplings are determined by 
fitting the t~ — > ii-k^-k°v t branching ratio and decay spectrum to recent 
data. Predictions for the partner decay r _ — > rj'-K^-K°v T and the low- 
energy behavior of the cross section a(e + e~ — > t]tt + tv~) are also discussed. 
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1 Introduction 



Tau decays represent an ideal benchmark where to analyze diverse topics in el- 
ementary particle physics [T|. In particular, semileptonic decay channels t~ — > 
H~ v T , where H is some hadronic state, allow a rather clean theoretical analysis 
of the hadronization of the V — A currents in presence of QCD interactions, 
since there is no hadron pollution to the leptonic current. Thus, these pro- 
cesses provide a suitable tool to find out intrinsic properties of the involved 
hadron resonances El HI El E] - In this work we concentrate on the analysis 
of r~ — > rj('\~Tr v T decays. For these channels the contributions of scalar and 
pseudoscalar resonances are expected to be negligible, since they turn out to be 
forbidden at tree level by symmetry arguments, such as G-parity conservation. 
In the limit of isospin symmetry the corresponding amplitudes arc driven by 
the vector current, allowing a precise study of the couplings in the odd- intrinsic 
parity sector. 

Concerning the theoretical description, it is well known that in the very low- 
energy domain (E <C M p , where M p is the mass of the p(770) meson) Chiral 
Perturbation Theory (xPT) [7] is the adequate tool to describe hadronic r de- 
cays [SJ. However, this approach fails when the invariant mass of the hadronic 
state becomes comparable with the mass of the lightest vector and/or axial- 
vector resonances, therefore a new strategy is needed in order to enlarge the 
domain of applicability of xPT to higher energies. One way out in this sense 
is to abandon the Lagrangian approach: one can model r decay amplitudes by 
taking the lowest order (LO) xPT results to fix the normalization of the form fac- 
tors at low energies, incorporating the dominant vector and axial-vector meson 
resonance exchanges by modulating the amplitudes with ad hoc Breit-Wigner 
functions [21 [31 HI El E] ■ However, it can be seen that in the low-energy limit this 
approach is in general not consistent with next-to- leading order (NLO) xPT [7], 
hence the usage of this procedure to reproduce QCD-ruled amplitudes is ques- 
tionable [TU1 QT] . An alternative approach is to include the lightest resonances 
as active degrees of freedom in the theory. This is done in the so-called Reso- 
nance Chiral Theory (RxT) [121 H31 [HI US] > where the inclusion of resonances 
is carried out together with an expansion in the inverse of the number of colors 
(Nc) [HI [17l [18l [19]. In fact, at the lowest order in the 1/JVc expansion one 
gets from QCD an effective theory that includes a spectrum of infinite zero- 
width states. The goal of RxT is to go beyond this limit, incorporating — in a 
way consistent with QCD symmetry requirements — a few finite- width resonant 
states that dominate the processes under study. 

A basic assumption in RxT is that the lightest resonant states are the dom- 
inant ones in low-energy phenomenology. In this way, for a given process it 
should be sufficient to introduce only the lightest resonance multiplet carrying 
the appropriate quantum numbers, while the inclusion of higher states can be 
carried out as a correction [201 HI]- On the other hand, the theory is built upon 
some fundamental QCD-based features: the effective interactions have to satisfy 
QCD symmetries, the low-energy behavior has to be consistent with xPT, and 
the asymptotic behavior of Green functions and associated form factors has to 
satisfy QCD constraints. These requirements imply several relations among the 
effective couplings that render the theory predictive. The aim of this work is to 
study within this framework the decays r _ — > rpt~ ir° v T , t~ — > r/' ir~~ ir a b> T , and 
the low-energy limit of the cross section <r(e + e~ — > r]ir + Tr~). 
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The article is organized as follows: In Sect. [5] we recall the fundamentals of 
Resonance Chiral Theory (see e.g. Ref. [52]) to explain our methodology. The 
relevant hadronic form factors for the decays under study are given in Sect. |5] 
In Sect. |H we derive the QCD-ruled high-energy constraints on the couplings, 
which reduce the number of unknowns to only four. In Sect. [5] we show that 
two of these unknowns can be bounded from other phcnomcnological studies 
performed within the same framework. In this way we end up with two unknown 
couplings, which appear to be highly correlated [231 l2~i] . The possible values 
of these couplings are analyzed by fitting experimental data on the differential 
decay distribution of t~ — » 7]Tr~ir°i' T and taking into account the present upper 
limit on the branching ratio for t~ — > 7]'ir~ir°v T . The low-energy behavior of 
the cross section cr(e + e~ — > rjTr + n~) is also discussed. Our conclusions are 
presented in Sect. [5] Finally, in Appendices A and B we analyze other possible 
contributions to the decay amplitudes and quote some useful isospin relations. 



2 Theoretical framework 

The effective Lagrangian in R^T is basically ruled by the approximate chiral 
symmetry of light-flavored QCD — which drives the interaction of light pseu- 
doscalar mesons — and the SU(3)v assignments of resonance multiplets [T2l[T4] . 
As we will see, for the processes under consideration it is possible to achieve a 
good agreement with present experimental data without the inclusion of excited 
multiplets. Moreover, it is seen that vector meson dominance (VMD) turns out 
to be a good approximation j!2j , since spin-zero resonance contributions vanish 
at tree level in the very accurate isospin symmetry limit (see App. [A|. In the 
case of r decays, owing to the relatively large r mass it occurs that several res- 
onances reach their on-shell condition when the amplitudes are integrated over 
the full phase space. The corresponding pole singularities can be regularized 
by including finite (energy- dependent) resonance widths, thus departing from 
the lowest order in the l/Nc expansion. Here we adopt the prescription in 
Ref. [25| . where energy- dependent resonance widths have been calculated in a 
well-defined way within RxT. 

We will work out r~ — > rf'' , n~'K (> i> T decays considering exact isospin symme- 
try. In this limit the processes are driven only by the vector current (see Sect. [3]), 
and appear to be dominated by the contributions of the p(770) resonance. The 
relevant effective Lagrangian in RxT reads : 

F 2 Fy 
£r x T = CwZW + £km + — (u^ + X +) + ^J^y^f+) 

r 7 4 5 

+ i^(V^u^) + E ^Mvjp + E d *°vvp + E ^^vppp . (!) 

^ i—1 ^ i—1 i=l ^ 

where all coupling constants are real, F and My being the pion decay constant 
and the mass of the lightest vector meson resonances, respectively. We follow 
here the notation in Rcfs. [TT1 [T2l l26f*1 . Accordingly, () stands for trace in flavor 

Ref. 1241 . two additional operators (Oyjp and Oyyp) have been found when the singlet 
(VVP) Green function is considered in addition to the octet one in the p 2 ~ m q ~ 1/Nq 
counting. In App. we show that they do not contribute to the hadronic tau decays studied 
here. 
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space, and u M , x+ an d /+ are defined by 



X± = x u ^ i u X u j 



(2) 



where u (U 



X and F£ r are 3x3 matrices that contain light pscudoscalar 



fields, current quark masses and external left and right currents, respectively. 
The matrix V^ v includes the lightest vector meson multiplet, and £jf in stands 
for the resonance kinetic term. The first term in Eq. ^ is the Wess-Zumino- 
Witten interaction Lagrangian [57J |2S], which governs the decay amplitudes 
studied here in the limit of low hadron momenta. The part of this interaction 
that contributes to the processes considered here reads 



C 



wzw 



iN c 

(L O R)) , 



s£ srr p 



(3) 



where E LjJ? are given by = Wd^U, = Ud^W, and Z Q and r Q are left 
and right external currents. Finally, the operators Cy JP , Oyvp an d Cvppp m 
Eq. (TTJ) are given by : 

VJP terms 



O 1 

^VJP 



^VJP 



^VJP 



FFP terms 



VJP 



Cvjp 

°VjP 
°VjP 



(4) 



^VVP 



O 2 



VVP 

3 

VVP 



CO 4 

^VVP 



(5) 
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VPPP terms 



Ovppp = e^ap ({V» u , u a u?} x-) , 

Ovppp = w(« a ^^X-)> (6) 

where = V M it„ + V^it^,. The covariant derivative V M involves pseudoscalar 
meson helds and P, r Q external currents. Its explicit expression can be found 
in Ref. [12]. 

The nonet of vector resonances V is described here using the antisymmetric 
tensor formulation. In the context of VMD [T3] , this is shown to be consistent 
with the usage of the xPT Lagrangian for light pseudoscalar mesons up to 0(p 2 ) 
in the even-intrinsic parity sector and up to 0(p 4 ) in the odd- intrinsic parity 
sector H5l. 



3 Form factors in r — > rj^ir tPv t 

In the Standard Model, r~ — ^ rf'' , K~'K v T decay amplitudes can be written as 

M = - ^| v ud u^Y 1 (i - 75K , (7) 

where V u d — cos Be is the relevant Cabibbo-Kobayashi-Maskawa mixing and 
Hfi is the hadron matrix element of the left-handed QCD current (V — A)^. In 
general, for a decay of a r lepton into three pseudoscalar mesons the hadronic 
tensor can be written as [5] 

(h 1 ( Pl )h 2 ( P2 )h 3 (p 3 )\(V - Ar\0) = F 1 A (Q 2 , Sl , S2 )V( l 

+ F 2 A (Q\ Sl , s 2 )V 2 u + iF 3 v (Q 2 , Sl , s 2 )Vg + F A (Q 2 , s 1} s 2 )Q», (8) 

where 

Vf = (<T - ^) (pi - p 3 )„ , Vf = (> ~ ^) (pa - pa), , 

= e^>i Q p 2/3 p 37 , <3 M = (Pi +P2+ PsT , Si = (Q - Pl f . 

(9) 

The upper indices in the form factors indicate the participating currents, either 
the axial- vector (A), or the vector one (V). The form factors F A and F A drive 
a transition to hadronic states with quantum numbers J p = 1 + , while F 3 V and 
F A correspond to outgoing states with J p = 1~ and J p = 0~, respectively. 
Let us focus on the amplitude for the transition r~ — > r?s(Pi) 7r_ (P2) 7r°(P3) v t , 
considering the limit of exact isospin symmetry. First of all, it is easy to see that 
for this process the axial- vector form factors vanish from G-parity conservation, 
therefore the dynamics will be essentially determined by the form factor F 3 V . 
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Figure 1: Topologies contributing to the final hadron state in r~ — > rf'' ir~ tt° v t 
decays in the Nc —> oo limit. Crossed circles indicate QCD vector current 
insertions. Single lines represent pscudoscalar mesons (77, tt) while double lines 
stand for p-rcsonance intermediate states. 



From the effective Lagrangian in Eq. (Q]), the diagrams that contribute to 
F§ are those represented in Fig. f, where single solid lines correspond to n 
and 77 mesons and double lines to the p(770) resonance. The corresponding 
contributions to the vector form factor read 

F V(a) N C , , 

3 {vs7T7r) 6\/6 7r 2 F 3 ' 

j-,V(b) 8Gy f r 2 

3 ^ ~ V3fzm v ^~t; [ci25Q - Ci25eSi 

+ 01235^ + 8c 3 {ml - m*)] , (11) 

E,V(c) 16 F V 1 , 2 2 2 \ 

- (<?! - ft,) 2mJ + g i5 ml] , (12) 

_ 8\/2 FyGy 1 1_ 

3 ^3 F 3 M 2 p -Q 2 M 2 p - Sl 

x [da (<3 2 + *i) + (du - da) m\ + 8d 2 (m* - m%)\ , (13) 
where we have defined 

C125 = ci - c 2 + c 5 , 
C1256 = ci - c 2 - c 5 + 2c 6 , 
C1235 = ci + c 2 + 8c 3 - c 5 , 

3123 = .9i + 2.92 - 93 , 
.945 = 2y 4 + 35 , 

di2 = di+8d 2 . (14) 

The amplitude for the r decay into the rjoir~ir° hadronic state can be read 
from Eqs. (TO]) to (HD by simply multiplying F^ (a ' b,c ' d) by y/2. Then, 
the matrix elements for the decays into the physical hadronic states r)ir~ir° and 
jfir~ir° can be obtained by considering rjs — ?7o mixing. Here we will consider 
a double angle mixing scheme |29j . which is consistent with the large- Nq ex- 
pansion [50]. Using a notation similar to that in Ref. [3T], the SU(3) octet and 
singlet fields are collected in a doublet rjg = (778,770), while the physical fields 
are included in r]p = (77,77'). These doublets are related by the transformation 
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r]B = (M) T r) P , where [3T] 

/ cos6» P (l - 5 s /2) + sm8 P 5 80 /2 - sin6» P (l - 5 /2) - cos<9 P S 80 /2 \ 

M = 

\ sin6» P (l - S 8 /2) - cos6 P S 80 /2 cos6» P (l - S /2) - sin6» P S 80 /2 J 

(15) 

In the framework of R^T, the parameters S 8 , So and S 8 q can be in fact derived 
from an effective Lagrangian that involves scalar resonances [13] . If the latter 
are organized in a U(3) matrix 5*, from the lowest order Lagrangian 



C s = c d (Su^) + c m (S X +) (16) 



one gets 



_ 8c d c m Mf 8c d c m Ml 8c d c m M| 

ds -^IfT^> d °~^ifT^> dm ~^if^> 



wherd^l 



M| = i(4M|-M 2 ) , 
M 2 = i(2M|- + A4 2 ) , 

M 8 2 o = -^(M^-M 2 ). (18) 

Here we take for M„ and Mr- the isospin averaged values of the pion and kaon 
masses, neglecting higher order corrections in the combined chiral and 1/Nc 
expansion. In addition we assume c d c m = F 2 /4 [521 EH, which is required by 
high-energy QCD in the Nq —> oo limit. Finally, from r\ — rf phenomenology 
we take M s ~ 0.980 GeV and 6 P = (-13.3 ± 0.5)° [35]. 

Given the form factors, F^(Q 2 , s\, S2), the spectral functions for the decays 
r~ — > ri^''K~'K Q v T are finally given by 

dr G%\V ud \ 2 (M 2 \ 2 1 / Q 2 \ 



dQ 2 128(2tt) 5 M t V Q 2 ) 3 V / 

/•(Vo 2 -™-) 2 /■t+(Q 2 ,s2) 
x / ds 2 / d*i W^B(g 2 ,si,s 2 ), (19) 

J(m„+m,) 2 Jt-(Q 2 ,s 2 ) 

where the relevant structure function W p [E] is defined by Wb(Q 2 , sx, S2) = 
^3 2 l-^TlQ 2 ' s ii s 2)| 2 and the limits of the integral over sx are 

i±(Q 2 ,.s 2 ) = ^ I (Q 2 +m 2 -2m 2 ) 2 

- A 1 /2(Q2 iS25TO 2 )TA l/ 2(TO 2 im 2 5S2) ] 2 | 5 (20) 

with A(o, b, c) = (a + b — c) 2 — Aab. We have neglected here the neutrino mass. 

2 The fully dominant contribution to the 7/ mass is not due to current quark masses but to 
the U{1)a anomaly 1321 . through the topological susceptibility of gluondynamics. Hence we 
keep Bp as a free parameter, to be fitted from phenomenology. 
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4 Short distance constraints on the couplings 



The above form factors depend on several combinations of coupling constants, 
besides the p mass and the pion decay constant. The values of these parame- 
ters are not provided by the effective theory, and their determination from the 
underlying QCD theory is still an open problem. However, one can get informa- 
tion on the effective couplings by assuming that the resonance region provides 
a bridge between the chiral and perturbative regimes, even when one does not 
include the full resonance spectrum [14] . This is implemented by matching 
the high-energy behavior of Green functions (or related form factors) evaluated 
within the resonance theory with asymptotic results obtained in perturbative 

qcd [mmnijuTjiMiiisisi. In p articular ' [t has been sh ° wn that the 

analysis of the two-point Green functions Ha, v [2] and the three-point Green 
function VVP of QCD currents (with the inclusion of only one multiplet of 
vector resonances) [35] leads to the following constraints in the Nc — > oo limit: 

i) By demanding that the two-pion vector form factor vanishes at high mo- 
mentum transfer one obtains the condition Fy Gy = F 2 |14| . 

ii) The analysis of the VVP Green function [26] leads to the following results 
for the couplings in Eqs. (HI]), (JT2J) and ([13]) : 



C125 = . 
C1235 = . 

C1256 = 



N C My 



327T 2 y/2Fy ' 

N c My F 2 



64tt 2 F 2 4F 2 
N c M v F 2 

V V 



d3 = -^-Fi + m- (21) 



On the other hand, it is possible to find additional constraints by requir- 
ing that the contributions of any intermediate hadronic state to the spectral 
function lmIly(Q 2 ) vanish in the limit Q 2 — > oo. This is a reasonable assump- 
tion, since from perturbative QCD ImIIy(Q 2 ) has to go to a constant value for 
Q 2 — > oo [37j , and the imaginary part of the two-point Green function can be 
understood as the sum of infinite intermediate hadronic states. Considering the 
intermediate rf^im hadronic states one gets the following constraints on the 
coupling constants: 



C125 = , 

C1256 = — 

d 3 = 

.9123 = , 



N C MyFy 



96tt 2 ^/2F 2 ' 

N c M v 
192tt 2 F 2 ' 



Nc My 

92 = I92^7^V- ^ 

It is worth to notice that relations (|22[) arc in agreement with those found in a 
similar analysis carried out for t decays into 2Kirv T [35] [3J5] and P~^v T (P 
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7r, K) states [3D]- Comparing with Eqs. (f^Tj) . we agree in the vanishing of C125, 
while the constraints for C1256 and d% cannot be simultaneously satisfied keeping 
agreement with their values in Eq. (J22). Moreover, as stated in Ref. [35], it is 
seen that the expected vanishing of the ^7*7 form factor at high-g 2 is obtained 
from Eqs. (|22p but not from Eqs. (|2~Tj) . In any case, numerically the differences 
are small, and the impact of these couplings on the observables is rather milojf|. 
Thus we choose to stick to our set of relations (|2"2"|) . using Eqs. (f2"T|) to fix the 
combinations C1235 and di2, not obtained within our study. In this way, the 
analysis of short distance constraints allows to reduce significantly the number 
of unknown coupling constants in the form factors F^ a ^ quoted in Eqs. (|10M13I) . 
To calculate the decay amplitudes for the processes t~ — > rf'> n~ tt° v t we end 
up with just four unknown parameters, namely Fy, C3, 545 and d^. As in the 
above mentioned analysis, we will take My = M p . 

5 Phenomenological analysis 

In order to carry out a phenomenological analysis of r~ — > rj^''n~ tPv t decays we 
take into account the available experimental information. In the case of r~ — > 
t[tt~ v r this includes the measured branching fraction BR(t" — > r\n~iJpv T ) = 
(1.39 ± 0.10) x 10 -3 [41], as well as the data on the corresponding spectral 
function obtained by Belle [Hj. The process r~ — > ri'-K~-K a v T has not been 
observed yet, hence we consider only the upper bound given by the PDG |41) . 
namely BR(t" — » r]''K~'K Q v T ) < 8.0 x 10~ 5 at 90% confidence level. 

As stated, the number of unknown parameters entering the vector form 
factor F% in R%T can be reduced to four by means of the short-distance con- 
straints obtained in Sect. 2] In addition, the values of Fy and 345 can be 
estimated within R^T from the phenomenological analysis of r~ — > (ttttt:)^ v t 
and to — > ir + 7T~ir°, respectively: the best fit to the — > (irinr)~ v T spectral 
function measured by ALEPH g3] corresponds to Fy = 0.180 GeV [II] with 
an estimated error of ~ 15?{|f|, while from the uj — > ir + ir~ir° branching ratio one 
gets (745 = —0.60 ± 0.02 [38] . In this way we are left with only two unknowns, 
namely the coupling constants C3 and c?2 ■ Our goal is to be able to describe the 
available experimental information just by fitting these two parameters. 

In Figure [5] we show the C3 — c?2 parameter region compatible with the PDG 
branching ratio for the mode t~ — > ?77r _ 7r°^ T at the level of one sigma. A 
large correlation between both couplings can be appreciated, in agreement with 
Refs. [531 124pl . Then, taking into account this allowed region for C3 and c^, 
we have carried out a fit to Belle data [55] for the t~ — > rjTr~TT°is T spectral 

3 The introduction of additional resonances has a different effect on the short-distance 
relations obtained from the VVP Green function and from the imaginary part of the vector- 
vector correlator. While all new contributions to the correlator are positive definite, this is not 
true for the VVP Green function, where cancellations are allowed. Thus the outcome of both 
procedures may be different when the spectrum is restricted to the lowest-lying resonances. 
A convergence of both results should be recovered if the full tower of excited resonances is 
taken into account. 

4 Some theoretical analyses lead to the value F v = V3F ~ 0.160 GeV [4011441145114611471 
[48] [49]. We have checked that a change of Fy within the range [0.160,0.180] GeV does not 
affect significantly the results presented in this section. 

5 fn particular, as noticed in Ref. | 24[ . there is an anticorrelation between C3 — in r~ — > 
rf n~ ir°v T and in associated radiative decays. Therefore, the combined study could improve 
the determination of these couplings. 
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function. We find two \ 2 minima, located at (03,^2) = (—0.018,0.45) and 
(c 3 ,d 2 ) = (0.035,-0.70), with x 2 / do1 = 3 - 80 and 4 - 34 , respectively, where the 
statistical error is about 10%. These values are indicated in Fig. [21 where black 
strips correspond to the C3 — d 2 regions that keep x 2 /dof within one unit far 
from the minima. The corresponding theoretical curves for the spectral function, 
together with experimental data, arc shown in Fig. [3J We have also carried out 
a fit to the normalized spectral function, obtaining that the preferred values 
for (c 3 ,d 2 ) remain almost unchanged, while x 2 /dof values get reduced to 3.0 
for both minima. Finally, in order to account for the theoretical error of the 
high-energy predictions for the couplings C1256, ^3, 32 and &\<z, we have also 
fitted the data allowing these coupling combinations to vary within ±1/3 of the 
values obtained from Eqs. (f2~Tj) and (|22{F1 . Noteworthy, the x 2 value does not get 
reduced, which can be taken as an indication that our short-distance relations 
(obtained at leading order in 1/Nq) lead to an appropriate effective Lagrangian 
to reproduce the experimental observations. 




-0.10 -0.05 0.00 0.05 0.10 



C 3 

Figure 2: Contour in the C3 — d,2 plane compatible with the branching ratio of the 
decay r~ — s- tyk~-k°v t at the level of one standard deviation. The black strips highlight 
the two regions that yield lowest values of \ 2 /dof (within one unit) from a fit to Belle 
data [42] on the corresponding spectral function. 



6 We have also considered nonvanishing values for the coupling C1235, which should be zero 
according to Eq. ill' I ft . Notwithstanding, have kept C125 and #123 equal to zero. Indeed, if 
C 125 7^ the Brodsky-Lepage behavior 1501 of the form factor is violated, and ImEIy goes, 
asymptotically, as Q 6 \og(Q 2 /My); if #123 7^ 0, the asymptotic growth goes as C(Q 6 ). Varying 
C1235 hi the range [—0.05, 0.05] does not improve the fit. 
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Considering the fitted values for C3 and d,2, we can calculate the correspond- 
ing predictions for the r~ — > rf ir~ 7r° 'v T branching ratio. The results are shown 
in Fig. |4j where we have taken C3 as the independent parameter. It is seen 
that the predictions are somewhat above the 90% confidence level upper bound 
quoted by the PDG, which is indicated by the shadowed region in the figure. 
However, the result corresponding to C3 = —0.018 turns out to be rather close 
to the upper bound; in fact, compatibility is achieved if the width of the C3 — c?2 
band is enlarged considering two standard deviations in the measured value of 
BR(r" — > r]ir~ ir° v T ) . Future, more precise measurements of the r~ — > i]'tt^tt°v t 
process should indicate whether our slight discrepancy arises from a weakness 
in the theoretical assumptions (e.g. treatment of r/g, — tjq mixing, effect of ex- 
cited resonances, SU(3) breaking terms in the Lagrangian [21]) or it just reflects 
an issue in the detection of this r decay mode. In this regard, we emphasize 
the importance of making global fits with unified and consistent treatments of 
all hadronic currents, in order to avoid cross-contamination between different 
hadronic tau decay channels from misunderstood backgrounds. The improve- 
ment in the most relevant hadronic matrix elements in TAUOLA [52] [S3] may 
be a key tool in this sense. 

Finally, our analysis can be used to predict a(e + e~ ?77r + 7r~) in the low 
energy region (conversely, one could in general use data on e + e _ annihilation 
into hadronic states to get predictions for the corresponding semileptonic tau 
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5.0x1 0" 4 




Figure 4: Prediction for the branching ratio of the decay r~ — > rj'-K~-K°v T consistent 
(within one sigma) with the r~ — > <ryK~i<Pv r branching ratio quoted by the PDG. The 
horizontal line, corresponding to BR(r~ — > -q'-K'-K^Ur) = 0.8 x 10 -4 , represents the 
current PDG bound. The notation for the black strip is the same as in Fig. [5] 



decays [5Tl l55fL The relation between this cross section and the • 
spectral function is detailed in Appendix [B] One gets 



T]1T 7T V T 



dT(r" -> ^-vrV) = 2 /(Q2)a(e+e - (23) 



dQ 2 

where f{Q 2 ) is given by 

/W J 384(27r)5M r ^ Q 2 ^ ^ + ^ M 2 )\4Sn) W ' 

In Fig. [5] we quote our predictions for the cross section, in comparison with 
low-energy e + e~ data obtained in various experiments Q We notice that al- 
though the ?/ meson decays to r]Tr + TT~ with a fraction of about 45%, there is no 
significant contamination from the chain er(e + e~ — > 7/7* ~ ^ ?/7r + 7r _ ) since, due 
to C parity, this occurs at NLO in powers of the electromagnetic coupling a. 
From the figure it is seen that our results are consistent with experimental data 
up to a center of mass energy of about 1.4 GeV. In fact, one should not expect 



7 A more elaborated dedicated approach, also based in RxT, has been developed for 
<r(e+e — —> ri/ir° 7r+7r — ) 1561 . 

8 We note that in this neutral current process there are additional contributions from new 
operators Oyjp and 0^ vp , see Ref. 1241 . which implies the introduction of two additional 
unknown couplings. However, these terms are suppressed in the large- Nq limit in the standard 
counting (see discussion in App.fSJ- 
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our treatment to be valid beyond this energy region, where effects of excited 
states should be sizeable and there is no phase space suppression as in r decay 
spectral functions. 



• DM1 
■ ND 

♦ DM2 

A CMD-2 

▼ BaBar 

- Our prediction 




1 1,2 1,4 



r)7m (GeV) 

Figure 5: Prediction for a(e + e~ — > T]n + %~) low-energy behavior from our analysis of 
r~ ?77r~7r i/ T decays, in comparison with DM1 [57], ND [58], DM2 [59], CMD-2 [60] 
and BaBar [61] data. 



6 Conclusions 

We have worked out the decays t~ — > iyk~tt°v t and t~ — > ii'ir~ir is T within 
the framework of Resonance Chiral Theory. The theoretical analysis has been 
based on the large-A'c expansion of QCD, the low-energy limit given by xPT 
and the appropriate asymptotic behavior of the form factors, which helps to 
fix most of the initially unknown effective couplings. Indeed, after taking into 
account information acquired in the previous related studies, t~ — > rj^ tt~ tt° v T 
amplitudes can be written in terms of only two unknown parameters. 

We have carried out a phenomenological analysis taking into account the 
experimental data for the branching ratio and the spectrum of the decay t~ — > 
ri7r~iT is Tl as well as the present upper bound for the t~ — > r)'ir~Tr°i/ T branch- 
ing fraction. A fit to the data allows to determine two preferred sets of val- 
ues for the unknown parameters C3 and d,2 in the effective Lagrangian, namely 
(c3,d2) = (—0.018,0.45) and (0.035,-0.88), which lead to a reasonable overall 
description of the spectrum. The former set seems to be favored by the prcdic- 
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tions for the branching ratio BR(t~ — > rf 'tt~ tt° v t ) , although in both cases the 
theoretical values appear to be somewhat above the present experimental upper 
bound. Finally, using isospin symmetry, these results can be used to get a pre- 
diction for the low-energy behavior of the cross section a(e + e~ — > ?77r + 7r~). The 
results are in good agreement with the available experimental information, and 
the approach can be useful |56j for the implementation of the related hadronic 
current in the PHOKHARA [HI] Monte Carlo generator. 

Our present results should be regarded as a first step in the study of the t~ — ^ 
7/')7r~7r ( V T decays in the framework of R^T. In the light of higher statistics 
for the ?77r~7r ^ T mode, or the observation of the decay — > T]' ~k~ it® v T , this 
description could be improved by considering e.g. the exchange of excited vector 
resonances, SU(3) breaking terms in the Lagrangian or revising the r]% — i]q 
mixing scheme. 
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A Contribution of spin-zero resonances and sin- 
glet terms 

In this Appendix we analyze both the contribution of spin-zero (scalar and pseu- 
doscalar) resonances and SU(3) singlet couplings to r~ — > r]^' ir~ tt° v T decays. 

A.l Scalar and pseudoscalar resonance exchange 

Discrete symmetries of QCD constrain the possible couplings in the effective 
Lagrangian. One of these symmetries is G-parity, which is exact in the SU(2) 
symmetry limit. With our conventions, the corresponding quantum numbers 
are: Ga^ = — 1, G% = +1, G n = +1, G„w = — 1, G/ / ff = +1, G ao = — 1, 
thus the final state t/tt^tt has G = +. As stated in Sect. 3, since the axial- 
vector weak current has G = — , only the vector current can contribute to the 
t~ — > rj-'^-K^TpVr decay amplitudes in this limit. The intermediate states foir~, 



13 



crir~ are also forbidden by G-parity conservation, which only leaves the channels 
aQ~7r° — > i]ir~ir° and aQ7r~ — > rjir°ir~. However, the vector current is J p = 1~, 
while the intermediate states clqit have parity P = (— 1) J+1 . Therefore, one can 
conclude that both scalar and pseudoscalar resonance contributions are strongly 
suppressed at tree level and can be safely neglected. 

A. 2 Contribution of double-trace terms 

In Ref. [53] two additional operators have been found with respect to those in 
Ref. (25J. Although these operators involve two traces, hence they are suppressed 
in the standard counting (in powers of p 2 ~ m 2 . and 1/Nc), it is seen that they 
become leading when a simultaneous counting in all three expansion parameters 
is carried out [30] . The operators read 

8 VJP = -i^M v J^e^ly^f? )log(det«) 

0\ V p = -id 5 M 2 ^e^ p „{V^Vn^og(detu), (A.l) 

where the tildes stand for u and / matrices that include the singlet term (and 
would contribute to the processes considered here through the 77 and rj' compo- 
nents of the rjQ meson). Once again the contribution of these operators vanishes, 
since the second operator only contributes to neutral current processes, while 
the first one leads to the contraction of symmetric and antisymmetric tensors 
in the r decay amplitudes. 



B Isospin relations 

In this Appendix we provide a derivation of Eq. (1231 , which allows to relate the 
r~ — > i]ir~ir°v T differential decay rate and the a(e + e~ — > i]ir + ir~) cross section. 

We work in the limit of SU(2) isospin symmetry, and we neglect Z-cxchangc 
contributions to the hadronic e + e~ cross-section, which is a safe approximation 
in the considered energy range. Thus this process will be driven by the vector 
current, via photon exchange. One expects to get a relation between this cross 
section and the vector current contribution to the decay of a tau lepton into the 
corresponding hadronic state. 

Since both rj% and 770 states are SU(2) singlets, we can compute isospin 
relations between 770,8"""" channels just by taking into account the isospin of 
7T7T states. Let us denote by T_ , To- the amplitudes ( ?77r7r|(ir^,w|0 ) and by 
T + _,T_ + ,T o the amplitudes -^(777r7r|(?ir^u — dT^d)\0), where the subscripts 
correspond to pion electric charges, and r\ can be either 770 or 77s. We obtain the 
relations 

-j= (T_ + T -) = ~ (T+_ + T_+ - 2T 00 ) = , 
7=2 (T °- T - 0) = -T2 {T - + - T+ - ) ' 

V3(T+_+T_ + -T 00 )=0, (B.l) 
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which lead to 

Too = , T+_ = -T_+ = T _ = -T_ • (B.2) 

Now let us consider the electromagnetic current. One can decompose it into 
1 = and 1=1 pieces: 

r£ m = 1 (2«y*u - - 5y 8 ) = r^ 0) + rfo , (b.3) 

where 

r o)) = J + *f* - 2 ^s) , j = i (sy'u - • (b.4) 

One can relate the amplitudes ( T77T7r|r M |0 ) for charge and isospin (Ar) \r]Trir) 
states by 

A+_ = -^=A l + -±=A , A_+ = —j^M + , A 00 = -~^ A o • ( B -5) 

Moreover the vanishing of the amplitude A2 implies 

2A 00 + + A-+ = . (B.6) 
In this way one obtains the following relations: 

A+- + A ao = -^=A 1 , A-+ + A 00 = > A i = ~~ h ~^~~ ± ' 

Ao = + ^ + = A °° = -VaAoo = f (A + - + , (B.7) 

which lead to 

|A+-+A-+| 2 + |A + _-A_ + | 2 = 2 (|A+_| 2 + |A_+| 2 ) = 4|A 00 | 2 +2|A 1 | 2 , (B.8) 

\A 1 \ i = \A + -\ a + \A-+\ 2 -2\A O0 \ 2 . (B.9) 
Thus corresponding cross sections are related by 

a (e+e - — > t]tttt^ \i=i = er (e+e~ — > ?77r+7r~) + cr (e+e _ — > ?77r~7r+) 

- 2 x 2 cr (e + e _ -> 7 ? 7r°7r°) 
~ 2cr (e+e~ ->• r/Tr+Tr") , (B.10) 

where the additional factor 2 in the first equation arises from the presence of 
identical particles in the final state. In the last line we have neglected the cross 
section to the 777r°7r° state, since it turns out to vanish at the lowest order in the 
electromagnetic coupling a owing to C-parity conservation. For the isoscalar 
part, from Eq. (|B.7|) we find 

a (e+e - -> -qnir) \ I=0 = 6cr (e+e - -> ?y7r°7r°) . (B.ll) 

Finally, one has 

1 ( T o- _ T -0) = ^ r o- = | ""- Jd | ) = A . (B.12) 
v 2 v 2 
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Taking into account that the e + e cross-section into three hadrons is given by 

<T e+e -^ hlh2h3 (Q 2 ) = ^6 / ds d * l^l 2 {~V^Vf) . (B.13) 

The cross-sections for the different modes read (|A + _| 2 = |^4i| 2 /2 = |A_ | 2 /2) 
a (e+e- h- ^tt") = |^ / ds dt |T_ | 2 (V 3 ,V^*) , 

<r(eV->^) = ^^/ dsd 4 |Too|2( ^ 3 "* ) ' (B - 14) 

where the additional factor 1 /2 in the second equation arises from the presence 
of identical particles in the final state. Thus one finally obtains 

— ^ ~ = f(Q 2 )v(e + e -> 777T7r)| /=1 

= 2/(Q 2 ) [er(e + e~ -> 77^+77") - 2(r(e + e" -> r/7r°7r )] 
~ 2f(Q 2 )a(e + e- -> TyTT+Tr") , (B.15) 

where f(Q 2 ) is the kincmatical factor given in Eq. (]24[) . 
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